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Abstract 



We discuss the residual isocurvature perturbations, fully-correlated with the 
curvature perturbation, that are automatic in the curvaton scenario if curva- 
ton decay is sufficiently late. We contrast these residual isocurvature pertur- 
bations with the generally un-correlated 'intrinsic' isocurvature perturbation 
generated by an additional field such as the axion. We present a general 
formula for the residual isocurvature perturbations, referring only to the gen- 
eration of the relevant quantity (Cold Dark Matter, baryon number or lepton 
number) in an unperturbed universe. Specific formulas for the residual isocur- 
vature CDM perturbation are given, for most of the commonly-considered 
CDM candidates. 
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I. INTRODUCTION 



The inhomogeneity of the early Universe first becomes measurable a few Hubble times be- 
fore cosmological scales enter the horizon. At this stage, well after Big Bang Nucleosynthesis 
(BBN), we know that the Universe consists primarily of photons, neutrinos, baryonic matter 
and Cold Dark Matter (CDM) [1]. The perturbations in the densities of these components 
are conventionally specified by the perturbation in the total energy density, characterised 
by a spatial curvature perturbation, (, and by three isocurvature perturbations 

with % =CDM, B or u, giving the inhomogeneity in the relevant number density per photon 
[2]. As we shall see, ( and the <Sj are constant from BBN to the approach of horizon entry. 
They set the initial conditions for the subsequent evolution of cosmological perturbations, 
and we refer to them collectively as the 'Primordial Density Perturbation'. 

Observation is consistent with the hypothesis that the isocurvature perturbations are 
absent, while the curvature perturbation is Gaussian with an almost perfectly flat spectrum 
[3-5]. On the other hand, spectral tilt, non-Gaussianity and isocurvature perturbations may 
all be present at some level. 

The Primordial Density Perturbation presumably originates during inflation, from the 
vacuum fluctuation of one or more light scalar fields. According to the usual hypothesis [2] , 
it originates solely from the vacuum fluctuation of the inflaton field, defined as the one whose 
value determines the end of inflation. This 'inflaton scenario' leads to definite expectations 
about the nature of the primordial density perturbation. There are no isocurvature pertur- 
bations, and non-gaussianity will almost certainly exist [6] only at the 10~ 5 level through 
second-order perturbations. 1 On the other hand, significant spectral tilt is expected, and 
significant running of the spectral index is quite possible [8,2]. It is also possible that the 
curvature perturbation is accompanied by gravitational waves of sufficient amplitude to af- 
fect the Cosmic Microwave Background (CMB) anisotropy [8,2]. These expectations have 
achieved the status of a 'standard model' of the early Universe, which is now used by almost 
all groups who analyse CMB and galaxy distribution measurements. 

It has been pointed out recently [9] (see also [10,11]) that the inflaton scenario is not the 
only possible mechanism to generate large scale structure from inflation. The primordial 
density perturbation may instead originate from the vacuum fluctuation of some 'curvaton' 
field, different from the inflaton field. The curvaton energy density has an isocurvature per- 
turbation, and this generates the curvature perturbation when the curvaton density comes 
to be a significant fraction of the total before the curvaton finally decays. 2 This 'curvaton 



x Bigger non-gaussianity of the \ 2 type can be generated in a multi-field model, but only at the 
expense of fine-tuning [7]. 

2 Recently, a completely different scheme has been proposed [12-14], in which the responsible field 
acts by perturbing the inflaton decay rate without ever contributing significantly to the energy 
density. 
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scenario' leads to completely different expectations about the nature of the primordial den- 
sity perturbation. Gravitational waves [15] (see also [16]), and probably also the spectral tilt 
will be negligible, but non-Gaussianity and isocurvature perturbations are perfectly possible 
[17]. 

The curvaton scenario has received a lot of attention [15,17-43] because it opens up new 
possibilities both for model-building and for observation. The present paper focusses on the 
possible isocurvature perturbations, following [17]. Because they have the same source as 
the curvature perturbation, they are fully correlated 

5i(x) = Si C(x) , (2) 

with Si numbers that depend on the physical model. Following [17], we call this kind of 
isocurvature perturbation residual. A given will be zero if the relevant quantity (CDM, 
baryon number or lepton number) is created after the curvaton decays. Otherwise it can be 
calculated within a given early-Universe scenario, and it is generally expected to be of order 
unity. Our main purpose is to derive specific formulas for the CDM case, for the various 
possibilities regarding the nature of the CDM that are commonly envisaged. 

In either the inflaton or the curvaton scenario, one can suppose that the vacuum fluctua- 
tion of one or more additional fields is also relevant. Such additional fields may generate ad- 
ditional isocurvature perturbations, which will in general be uncorrelated with the curvature 
perturbation and at best only partially correlated. We shall call these intrinsic isocurvature 
perturbations. In general, intrinsic isocurvature perturbations and the curvature perturba- 
tion are determined by different sectors of the underlying particle theory. This means that in 
general, there is no reason to expect the magnitude of an intrinsic isocurvature perturbation 
to be comparable with that of the curvature perturbation. In a particular model though, it 
may turn out that a single sector of the theory accounts for both the curvature perturbation 
and an intrinsic isocurvature perturbation, so that their magnitudes are comparable. (See 
[36] for the only example known so far, arising within the curvaton scenario.) 

The plan of the paper is as follows. In Section II we recall the general description of the 
Primordial Density Perturbation, in terms of quantities which are conserved on super-horizon 
scales. In Section III we discuss how isocurvature perturbations can originate, first in the 
inflaton scenario and then in the curvaton scenario, giving a general formula for the residual 
isocurvature perturbations. In Section IV we apply our formula, to evaluate the residual 
CDM isocurvature perturbation for each of the commonly-envisaged CDM possibilities. We 
conclude in Section V. 

II. CONSERVED PERTURBATIONS 

The discussion of the evolution of perturbations while they are outside the horizon is fa- 
cilitated by the existence of perturbations that are, under suitable circumstances, conserved 
[44,42]. The conserved perturbations are conveniently defined with reference to the slices 
of spacetime which have zero intrinsic curvature perturbation (spatially flat slices), because 
on super-horizon scales the local expansion of the Universe between such slices is uniform 
[44,42]. By virtue of this uniformity, a perturbation will be conserved on super-horizon scales 
if it is constructed according to the following recipe [42]. Consider a quantity /(x, t), with t 
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the coordinate time labelling spatially flat slices. If / decreases (or increases) monotonically 
with the comoving volume V according to an equation of the form 



then the perturbation 



V§£ = </(/), (3) 



X f = ~H 6 -l (4) 



is conserved [42], where Sf is to be evaluated on spatially flat slices. In practice we find 
conserved quantities of three types. 

a. The total energy density perturbation. The quantity 

(e-^ = A (5) 
p p + P KJ 

is conserved provided that the pressure P of the Universe is a unique function of the energy 
density p, or equivalently if the perturbations satisfy the adiabatic condition 

5P/5p = P/p . (6) 

More general, if the pressure is non-adiabatic, 

C = --^p«W, (7) 
p + P 

where 

5P nad = 5P--5p. (8) 
P 

b. Separate energy density perturbations. Let pi and p refer to some component which 
by itself satisfies the adiabatic condition Eq. (6), and which does not exchange energy with 
any other component. Then the quantity 

Pi Pi + Pi 

is conserved. If the component is radiation (Pj = Pi/3), 

1 5 pi 

Q = - — (radiation) , (10) 

4 Pi 

and if it is matter (Pj = 0) 

d=\ 5 -^ (matter). (11) 
3 Pi 

The conservation of Q in these cases is an immediate consequence of the dependence of the 
local energy densities on comoving volume V (namely pi oc V -1 for matter and Q oc V~ 4//3 
for radiation). 
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c. Number density perturbation. Let Hi be any conserved number density. Then, since 
rii oc V -1 the quantity 

^l-tr (12) 

is conserved. 



A. The primordial density perturbation 

The epoch T ~ 10 MeV, marking the beginning of of the Big-Bang Nucleosynthesis 
process, is the earliest one of which we have definite knowledge. After T falls below 1 MeV, 
the energy density is dominated by radiation in the form of photons and neutrinos, and there 
is also baryonic matter and presumably Cold Dark Matter (CDM); 3 The number of particles 
is conserved for each component, and there is no exchange of energy between components. 
As a result the following quantities are conserved until the approach of horizon entry: 



_ l^Pcdm 
Scdm — „ 

" Pcdm 


(13) 


> !<W 
3 Pb 


(14) 


^ U ~ A 

4 p v 


(15) 


C7 " 4 p 7 


(16) 



The three isocurvature perturbations Eq. (1) are given by 

Si = 3(0 - C 7 ) , (17) 

and they too are conserved until the approach of horizon entry. The total energy den- 
sity/curvature perturbation, Eq. (5), is given by 

P 

^ (1 - fu)Cr + fuCu , (19) 

during the primordial era when the matter density is negligible and f u = p v j p 7 is a constant. 
The curvature perturbation ( is thus constant during the radiation dominated era to high 
accuracy until the approach of horizon entry. 

The smallest cosmological scale, enclosing say M ~ 1O 6 M , enters the horizon at T ~ 
1 keV. The curvature perturbation ( and the isocurvature perturbations St therefore have 



3 If non-relativistic particles decaying well after nucleosynthesis have a significant effect, we are 
assuming here that the perturbation in their energy density is adiabatic and hence need not be 
considered separately. 
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constant values on all cosmological scales in the regime 1 keV < T < 1 MeV. We are 
calling them collectively the Primordial Density Perturbation, and through the Einstein 
field equation and the Boltzmann hierarchy, they specify the subsequent linear evolution 
of the entire set of cosmological perturbations [2,45]. They are therefore determined by 
observation. 

It is found from observation [3-5] that the primordial curvature perturbation is almost 

i 

Gaussian, with an almost scale-independent spectrum given by V£ ~ 5 x 10 -5 . Current 
observations provide only upper bounds on the isocurvature perturbations [46,47,31,4,48,49], 
which depend on the degree of correlation between them and the curvature perturbation. 
In this paper we are concerned with the fully-correlated perturbations that can be produced 
only in the curvaton scenario. For them, the ratios «Sj(x)/£(x) are constants determined by 
the particle physics model. 

The observational bound in the case of sb = 5b( x )/C( x ) is [31] at 95% confidence level 

-0.53 < s B < 0.43. (20) 

The bound on s cdm = S , c d m ( x )/C( x ) is tighter by precisely the factor Q B /Q cdm [31], and 
taking that ratio as 1/6 we get 

-0.09 < Scdm < 0.07. (21) 

The observational bound on s u is [49] 

-0.14 < s v < 0.47. (22) 

This bound has been calculated for negligible lepton asymmetry, whereas according to 
Eq. (29) below a non-zero lepton asymmetry is required [17] to produce an neutrino isocur- 
vature density perturbation. Nevertheless, because the lepton asymmetry is small the bound 
may be expected to provide a good approximation. 



B. Conserved number densities 

To describe the origin of the isocurvature perturbations in the early Universe, we need 
to consider [17,42] the CDM number density n c d m , the density of baryon number n B and 
the density of lepton number n^. Their perturbations may be characterised by [17] 

C, = -a— , (23) 

Hi 

with % =cdm, B or L. We can take the epoch of creation of CDM, baryon number or lepton 
number to be the epoch after which the relevant quantity (CDM particle number, baryon 
number or lepton number) is conserved. Each quantity Q = Sni/3rii is then conserved after 
the relevant epoch of creation [42]. 

In the case of CDM, £ c d m in Eq. (13) and £cdm Eq. (23) will usually be equal from the 
moment of CDM creation, and in any case are equal by the primordial epoch. The exception 
is the case of axionic CDM though, where the axion mass is initially temperature-dependent 
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so that only £ c d m is conserved after creation. In that case, Ccdm, will be equal to Ccdm after 
the mass becomes constant. 

It follows from Eqs. (17) and (19) that 

<S c dm = 3 (Ccdm — C7) (24) 

= 3(Ccdm - C) • (25) 

The final equality assumes that S v = 0, but the modification in the case of nonzero S v is 
quite straightforward [17] and there is no need to consider it here. 

In the case of baryonic matter, there is no useful definition of (b before the quark-hadron 
transition, since baryon number is not carried by a particular energy density. Afterwards, 
Cb = Cb leading to 

Sb = 3(Cb - Q (26) 
= 3(Cb - C) • (27) 

In the case of neutrinos, ( u is of course always defined provided that neutrinos exist but 
it is a useful quantity only after neutrinos freeze out of thermal equilibrium. After that 
epoch [17], 



5„ = 
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(!) (a-c) (28) 

(!)VC)- (29) 

In these expressions, (l is the perturbed net lepton number, Siil/Sui, and £ is the neutrino 
asymmetry parameter, constrained by nucleosynthesis to |£] < 0.07. The second equality is 
a good approximation [17] by virtue of the bound on |£|. 

In the curvaton scenario, (l is negligible compared with ( if the neutrino asymmetry is 
created before the curvaton contributes significantly to the energy density [17]. Then, the 
present observational bound Eq. (22) leads to — £ < 0.27, which is competitive with the 
nucleosynthesis bound. 



III. INTRINSIC AND RESIDUAL ISOCURVATURE PERTURBATIONS 
A. The inflationary initial condition 

Inflation is supposed to set the initial condition for the subsequent evolution of the 
Universe, through the values of the light scalar fields which exist at the end of inflation. We 
are here defining a field as light if the second derivative of the potential in that direction is 
much less than H 2 . This means that a light field has a flat spectrum of perturbations with 
spectrum (H^/27r) 2 , with the star denoting the epoch of horizon exit. 

To handle the perturbed Universe, we adopt what has been called the separate universe 
picture [44,2], which means that the local evolution of regions within our presently observable 
Universe is the same as that of an unperturbed FRW universe. One assumes that the 
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Universe is smoothed on the shortest relevant scale. At least for cosmological scales the 
separate universe picture seems very likely to be valid, since such scales are far larger than 
the Hubble scale in the very early Universe. 4 As a result, physical processes in the very early 
Universe are unlikely to distinguish between the shortest cosmological scale (say, 1 Mpc) and 
the longest one corresponding to the present Hubble scale (~ 10 4 Mpc). And the separate 
universe assumption is certainly valid on the latter scale, or else the whole concept of an 
using an unperturbed FRW background to describe our observable Universe makes no sense. 

For an unperturbed FRW universe, a generic quantity g (such as energy density, pressure, 
number density or a scalar field) has an evolution of the form 

</(0i, 02,-.., , (30) 

where 0, are the values of the light fields specified at some initial epoch. We are specifying 
the epoch at which g is to be evaluated by the amount of expansion 

N=\ ln(V/V initial ) = ln(a/a initiaI ) = f Hdt . (31) 

<J -^initial 

In the separate universe picture, the locally-defined generic quantity g has the same func- 
tional form as in the unperturbed Universe, Eq. (30), but with position-dependent argu- 
ments; 

0(0 1 (x),0 2 (x),...,7V(x,t)), (32) 

where (x, t) are the coordinates (gauge) used to describe the perturbations, and N is the 
locally-defined defined quantity. 5 

The first-order perturbation can be written as 

^ = E^W + ^iV(x,t). (33) 

For a given wavenumber k/a one can choose the initial time to be a few Hubble times after 
horizon exit during inflation. Then the S(f>i(k) are uncorrelated gaussian quantities, with an 
almost flat spectrum 

V(k) = (Hj2n) 2 , (34) 

where the star denotes the epoch of horizon exit aH = k. 

Any conserved quantity Xf, defined in Eq. (4), can be calculated by evaluating 5f on 
a spatially flat time-slice as a function of the initial field perturbations . This is separated 



4 In [44,2] the separate universe picture at a given epoch is taken to be valid for any smoothing 
scale bigger than the Hubble distance. That is a stronger assumption, which might not be valid in 
the presence of short-distance phenomena like preheating [51]. 

5 The final expression of Eq. (31) is now / Hdr where H is the locally-defined quantity Hubble 
expansion and r(t) is the local proper time. 
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from an initial spatially flat time-slice during inflation by a uniform expansion, i.e., SN = 0, 
which yields from Eq. (33) 

*' = -7 v = -7© w * w ' (35) 

In particular we can write the total curvature perturbation as 

5<Mx) , (36) 

N 

This is a generalisation of the approach advocated by Sasaki and Stewart [52] as the 
simplest way to calculate the curvature perturbation in multi-field models. They calculated 
( = 5N in the particular case where N is the local expansion between an initial spatially 
flat time-slice and a final uniform-density slice [44]. Taking g = p and requiring 5p = in 
Eq. (33) reproduces this expression: 

C = W=-Eg) (37) 

where we have used dN / dp = H/ p along the background solution. 

The general expression Eq. (33) has not been written down before. But it is all that is 
required in the separate universe picture to calculate the conserved quantities Q Eq. (23) 
produced from field fluctuations during inflation, and hence describe the Primordial Density 
Perturbation. In what follows we will be interested in the case of isocurvature perturbations 
corresponding to conserved number densities n,, and our main result will be Eq. (46), which 
along with Eq. (66) is a further example of Eq. (33). 




B. The inflaton scenario 

During single-field slow-roll inflation, 6 the inflaton field satisfies 3/70 = — V ((/>), which 
can be integrated to give as a unique function <f>(N) of the integrated expansion N, up to 
a constant of integration. 

In the the inflaton scenario, <fi sets the initial condition for both the energy density and 
the pressure. As a result, the pressure perturbation is adiabatic, making ( constant. Using 
Eq. (5) with bp ~ V'8(f) one finds that a few Hubble-times after horizon exit the well-known 
formula 

C--f^, (38) 

where the right hand side can be evaluated at horizon exit. In the inflaton scenario, this 
value is maintained until the primordial epoch. 



6 We adopt the now-standard terminology, where 'single-field' means that there is an essentially 
unique slow-roll trajectory while 'multi-field' means that there is a family of slow-roll trajectories. 
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A primordial isocurvature perturbation Si can be generated in the inflaton scenario if 
the relevant number density rij at the time of its creation depends on some scalar field x as 
well as the inflaton field; 



We shall call this kind of isocurvature perturbation an intrinsic one, as opposed to the 
'residual' one that can appear only in the curvaton scenario. The known cases where it is 
possible are when the CDM comes initially from an oscillating field (axion or Wimpzilla), 
or when there is Affleck-Dine baryo- or lepto-genesis. The qualification 'possible' is here 
crucial, because the axion or Affleck-Dine field has to be light during inflation so as to 
acquire a perturbation, and also because that perturbation has to survive until n; is created. 
This does not happen if the CDM axions come mostly from the oscillation of Peccei-Quinn 
cosmic strings, or if the Affleck-Dine field has the generic effective mass-squared of order 



The possibilities for obtaining an intrinsic isocurvature perturbation are very limited, 
so that, for example, there is no known way of obtaining an intrinsic CDM isocurvature 
perturbation if the CDM consists of the relic abundance of the Lightest Supersymmetric 
Particle obtained when it decouples from thermal equilibrium. Even when there is an in- 
trinsic isocurvature perturbation, there is in general no reason to expect its magnitude is 
comparable with that of the adiabatic perturbation. This is because the magnitudes of the 
two perturbations depend on different and in general unrelated parameters, as a result of 
the fact that the perturbations are caused by different fields. 

A final and very important property of an intrinsic isocurvature perturbation is that it 
is uncorrelated with the adiabatic perturbation. Again, this comes from the fact that two 
different fields are involved. 

With one possible exception, the above discussion applies also to multi-field inflation, 
if we define the inflaton as the specific linear combination of fields corresponding to the 
inflaton trajectory at the end of inflation. The possible exception concerns the correlation; 
there is correlation in the exceptional case that the combination of fields orthogonal to the 
inflaton was responsible for the isocurvature perturbation [50]. However, no example of 
this has been proposed in the context of particle physics; in particular, no model has been 
exhibited where this combination is the axion or an Affleck-Dine field. 



In the curvaton scenario, the value Eq. (38) of ( generated by the inflaton is supposed 
to be negligible compared with the observed value of order 10~ 5 , and to remain negligible 
until the curvaton field a starts to oscillate. The density p a is supposed to be negligible 
when the oscillation starts, but during radiation domination it grows like a(t) relative to the 
radiation density Until the curvaton starts to decay, the energy densities of the radiation 
and the curvaton are decoupled, so that ( T is constant and negligible, while ( a is constant 
but not negligible. From Eq. (36) we then have 



(39) 



±H 2 . 



C. The curvaton scenario 




(40) 
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where ( a = —H5a/a and 




Pa 3p CT 



(41) 



p 4p r + 3p CT ' 



is a growing function of time. 

When curvaton decay is complete, the Universe is supposed to be thermalised, so that 
the temperature determines the entire future evolution of the energy density and pressure. 
(Non-thermalised components such as CDM are supposed to give a negligible contribution 
to these quantities until at least the primordial epoch.) As a result, the pressure in each 
separate universe is once again a practically unique function of energy density, leading to a 
practically constant value of ( given by Eq. (43). Following [17,28] we define 



where ( on the left-hand-side is evaluated after curvaton decay is complete, and ( a on the 
right-hand-side is evaluated before the decay starts. 
In the sudden decay approximation we have 



Hence we have r ~ (po-/p)dcc- Equation (43) gives a good approximation (to within 10%) 
to the precise numerical solution [28]. Note that r is found to be always less than (or equal 
to) one and, in particular, r — > 1 in the limit p a j p — > 1. 

In the curvaton scenario there might be an intrinsic primordial isocurvature perturba- 
tion due to some other field \i different from both the inflaton <fi and the curvaton a. If 
present, it will have the features discussed already in the case of isocurvature perturbations 
for the inflaton scenario. It will generally be uncorrelated with the residual isocurvature 
perturbation and accordingly we ignore it from now on. 

We are interested in a different type of primordial isocurvature perturbation, which we 
term residual [17]. A residual isocurvature perturbation can arise because in the curvaton 
case, the separate universes are not identical while the curvaton is oscillating. This is because 
they are characterised by the two independent local energy densities p^ and p a . Thus, a 
residual isocurvature perturbation Si will inevitably be present at some level unless the 
relevant number density is created after the curvaton has decayed. 

In the sudden-decay approximation, it is easy to arrive at a general formula for the 
residual isocurvature perturbations. Note first that before curvaton decay, p r is supposed 
to be uniform on the flat slicing. The perturbation of the number density nj(p r ,p CT ) on this 
slicing is therefore just given by 



r = C/C 



(42) 



(43) 




(44) 



As a result the perturbation, Eq, 



(12) 



can be written as 




(45) 
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Using Eq. (42) for the total curvature perturbation after curvaton decay, we have 



) 



(46) 



C r n \dp a 



This expression can be evaluated at any epoch after the conservation of rij but before curva- 
ton decay, since both and p a are proportional to a~ 3 . Courtesy of the separate universe 
assumption, its evaluation requires only a knowledge of the function rii(p T , p a ) in an un- 
perturbed universe. This is similar in spirit to the calculation in [52] of the curvature 
perturbation produced by multi-field inflation. The separate universe assumption can be 
verified by explicit calculation in all three of the CDM cases that we shall deal with in the 
next section, since the situation there is sufficiently simple, just as it can be verified in the 
case of multi-field inflation. 

In the case where the CDM or baryon number is created before the curvaton decays, the 
residual CDM or baryonic matter isocurvature perturbation is thus given by Eq. (2) with 



In the opposite case where the curvaton decays and its decay products are thermalised before 
the CDM or baryon number is created, the CDM or baryon number at creation is simply a 
function of the local energy density, as in the inflaton scenario, and we have Sj = 0. 

IV. THE RESIDUAL CDM ISOCURVATURE PERTURBATION 

For the rest of this paper, we focus on the residual CDM isocurvature perturbation, 
which is present in the curvaton scenario if the CDM is created before the curvaton has 
decayed. We shall work out Eq. (47) in the cases that the CDM (i) is a WIMP (ii) is an 
axion and (iii) comes from an oscillating field with constant mass, such as a Wimpzilla. We 
shall also compare our results with present observational bounds. 

Before looking at the individual cases though, we recall two results already considered 
in [17]. The first result holds in the case that the CDM is created at an epoch when (p a /p) 
is negligible compared with the value when the curvaton decays, which means in particular 
that the CDM must be created well before the curvaton decays. If the curvaton density has 
a negligible effect on n cdm , Eq. (47) becomes 



This is too big to be compatible with the current observational bound given in Eq. (21), 
and is independent of the nature of the CDM. It follows that in the curvaton scenario, CDM 
cannot be created significantly before the curvature perturbation achieves its final value. In 
particular, CDM cannot be created just after inflation, as would be the case if the CDM 
consisted of heavy weakly interacting particles created from the vacuum. (This kind of CDM 
was originally envisaged in [53,54], and later [55] called Wimpzilla CDM.) 




(48) 



(47) 



= -3. 



(49) 



12 



The other result holds when the CDM is created by the decay of the curvaton itself. The 
epoch of CDM creation then corresponds to the epoch when the curvaton decay is complete. 
The resulting local CDM density is a fixed multiple of the curvaton number density well 
before decay. The fractional perturbations are thus equal and hence 

Ccdm = C- (50) 

In terms of r, defined in Eq. (42), the prediction Eq. (50) leads to 




This is compatible with Eq. (21) only if 1 — r < 0.03, which means that the curvaton has 
to dominate before it decays. Note that r < 1 [28] and hence in this case there is a positive 
correlation between the primordial curvature perturbation ( and isocurvature perturbation 
S, i.e., s > 0. 



A. Weakly interacting massive particles 

If the CDM consists of weakly interacting massive particles (WIMPs), it is initially in 
thermal equilibrium. The number density of the WIMPs is conserved only after the epoch 
when the WIMPs fall out of thermal equilibrium (the epoch of freezout). In the present 
context that epoch should be taken as the one when the CDM is created. In contrast with 
the axion case that we shall discuss next, no scalar field can be involved in the creation of 
WIMP CDM and therefore there can be no intrinsic isocurvature perturbation. 

The calculation of the residual isocurvature perturbation after freeze-out in terms of 
gauge-invariant linear perturbations is given in an Appendix A. However we shall show here 
how the same result can be obtained working only in terms of the unperturbed background 
equations, via Eq. (47). 

In either case we work in the approximation of sudden freeze-out, denoting the moment 
of freeze-out by a star. This moment is determined by an equation of the form 

r(T*)/H* = K, (52) 

where K is of order 1 and T is the interaction rate per CDM particle. The temperature 
dependence of the interaction rate given in terms of a dimensionless parameter 

d\nT 
din 7" 

with a > 2 if there is initially thermal equilibrium at high temperatures. 
The number density of non-relativistic particles (m > T*) at freeze-out is 

n * =9 *{^t e ~ m/T% (54) 

where m is the mass and the internal degrees of freedom of the WIMP. The annihilation 
cross-section for non-relativistic WIMPs is weakly dependent on temperature, with {(Tav) oc 
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(53) 



T p l 2 where p = for s-wave annihilation and p = 2 for p-wave annihilation [1] . The rate of 
change of the interaction rate T is thus mainly determined by the rapidly decreasing number 
density if m > T, yielding 



T = n(a A v)ocT 3 / 2+ ^e- m / T , 



(55) 



and hence a = (m/T) + (3/2) + (p/2). 

After freeze-out interactions are negligible and the number density is diluted as the 
universe expands, so we have n = n*(a*/a) 3 . During any era when the number of effective 
species is constant, we have a oc 1/T and hence 



n = ?% 



(T/T*) 3 = ^(m/27rT,)le- m / T *T 3 , 



(56) 



This is the expression that we want to insert into Eq. (47). The partial derivative is at 
fixed p r and therefore at fixed temperature, T < T*, so that Eq. (47) can be written 



, o _ 3 (pa dn 

■Scdm + <J — ,-, 

r \n Op a 



(57) 



The local CDM number density after decoupling varies with the local curvaton density 
because the curvaton density at freeze-out, p aif , affects the expansion rate, in Eq. (52), 
and hence the local freeze-out temperature, Hence we can write 



' p a dn \ 
n dp a J T 



m 



' Pa dTA 
T* dpa) T 



(58) 



where we have used Eq. (56) to obtain (dn/dT*) T = {m/T 2 — 3/2T*)n. 

Assuming that the curvaton is oscillating at freeze-out (i.e., the curvaton mass is greater 
than the Hubble rate, m a > H*) we have 



Pa Pa* 



Po 



rp \ 3 

tZ 



(59) 



and hence 



Pa dTA 

T* 9 Pa J rp 



Pa dT* 



Finally then we must determine dp a */dT* from Eq. (52), which gives 



dK 



dT* 

~f7 



1 dp* 



(a* - 2(1 - 



r(r.) 



- -n 



dp 



2 Pa 



= 0. 



(60) 



(61) 



and hence 
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% dp a , _ 2(t*« -2) + 4fl c 
Pa dT* 



(62) 



Combining Eqs. (57), (58), (60) and (62) we then obtain the residual isocurvature per- 
turbation, 



^cdm 



m 



27 2(a* 



2) + n 



- i 



(63) 



In the usual case that freeze-out occurs during radiation domination (requiring <C 1) 
with the WIMP a neutralino, then m/T* is of order 20. In this case the annihilation cross- 
section is approximately constant while the number density drops rapidly and hence 



and Eq. (63) reduces to 



' d In n " 
dlnT 



m 



(64) 



^cdm 



2r 



- 1 



(65) 



The first term in the bracket is thus roughly of order /Q a dec (using the sudden-decay 
approximation Eq. (43) which gives r ~ fio-dec)- If freezout occurs well before the curvaton 
decays and while the curvaton accounts for only a small fraction of the energy density (i.e., 
for <C Q adcc ), this makes s c d m — —3 in accordance with the general expectation [17] for 
this regime, which violates the observational bound. 

On the other hand, if the curvaton dominates the energy density before freeze-out so that 
f2 CT * ~ ^o-dcc ~ 1 (or if freeze-out occurs close to curvaton decay with fl ~ Q(j dec ^ 1); 
there could be a regime where s c d m is small enough to satisfy observational bounds. Note 
that if the CDM freezes-out when the curvaton has come to dominate the energy density 
(Q a * ~ 1) but before it decays, then the CDM abundance relative to the radiation density 
will be significantly diluted when the curvaton decays. Thus m/T* could be significantly 
smaller than usually assumed. In order to determine the allowed regime more precisely 
would require more detailed numerical modeling going beyond the sudden freeze-out and 
sudden decay approximations that we have used here. 



B. Axion CDM 

If the CDM consists of axions, it is created when the rising axion mass m(T) becomes 
equal to the falling Hubble parameter H because only then does the approximately homoge- 
neous axion field become free to oscillate. If axionic strings are present, the axions radiated 
from them will probably be the dominant creation mechanism. In that case there can be 
no intrinsic CDM isocurvature perturbation. If instead such strings are absent, the axion 
field x(x, t) oscillates independently at each comoving position, with an initial value that 
is the same as it was during inflation. (This is the separate universe assumption, which in 
the present case can be verified by explicit calculation [56] from the axion field equation.) 
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Then, the axion density is proportional to the square of the initial axion field value, and a 
perturbation in this value will generate an intrinsic isocurvature perturbation 



'cdm 



2Sx 

X 



(66) 



We are here concerned instead with the residual isocurvature perturbation, which is 
always present if the curvaton decays after the epoch of axion creation. We work in the 
approximation that the oscillation starts suddenly and is immediately harmonic, which is 
known to be reasonable in most of parameter space [56]. Then, if axion creation occurs 
during radiation domination the calculation is very similar to the one that we gave for the 
WIMP case. The freezout epoch is determined by 

m(T*)/H* = K, (67) 

with K conventionally set equal to 1 exactly. The increasing mass is given accurately by 

mocT^, (68) 

with (5 = —3.7 [57]. Assuming that there is no intrinsic perturbation, the axion number 
density n* at creation is proportional to m(T*). Repeating the calculation of the WIMP 
case we find for the axion case 



Scdm 



a 



3 + p 



tt adec 2(3 + A -n c 



(69) 



In contrast with Eq. (63), this always gives s cdm ~ —3. We conclude that if the CDM 
consists of axions, the curvaton cannot decay after the CDM is created at T ~ GeV. 



C. CDM from an oscillating field of fixed mass 

The previous calculation does not apply if axion creation occurs at low temperatures, 
and in particular if it occurs while the curvaton density dominates. In such a case, the 
temperature of the sub-dominant radiation will have a negligible effect on the axion mass, 
which will therefore have its vacuum value. 

More generally, the CDM might consist of particles corresponding to the oscillation 
of some scalar field other than the axion, in which case the mass will generally be at the 
vacuum value whether or not creation takes place during matter domination. This possibility 
is envisaged in the rather attractive proposal of Moroi and Randall [58] , whereby a modulus 
decays before nucleosynthesis into both baryonic matter and CDM. The decay occurs because 
the modulus mass is 10 to 100 TeV instead of the usual 1 TeV or so. In contrast with earlier 
proposals though, the high mass is explained naturally, arising from the fact that in this 
scenario SUSY-breaking is anomaly-mediated as opposed to gravity-mediated. The anomaly- 
mediation also ensures that matter will not be over-produced by the decay. Of course, this 
idea has been worked out only under within the inflaton scenario and it remains to be seen 
whether it can be viable in the curvaton scenario, with the oscillation starting before the 
curvaton decays as we are envisaging. 
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Without focusing on any particular scenario, we consider the residual CDM isocurvature 
perturbation that is produced if the CDM does originate from the oscillation of some field 
with fixed mass. With constant mass, corresponding to (5 — 0, Eq. (69) can be written 

- l) , (70) 




^ct ,dec 

where / is defined in Eq. (41). Using Eq. (40) this equation corresponds actually to 

Ccdm = C* • (71) 

This in turn is simply a consequence of the fact that the oscillation starts on a slice of 
uniform energy density, corresponding to local Hubble parameter H = m. 

If the oscillation starts well before the curvaton decays, and while the curvaton accounts 
for only a small fraction of the energy density, this again gives s c d m — — 3 in contradiction 
with observation. However, in contrast with the cases of WIMP and axionic CDM, s cdm in 
the present case vanishes in the limit that the curvaton completely dominates the energy 
density when the CDM is created (ie. when the CDM oscillation starts). Using Eq. (70) one 
finds in this regime 

Srf „ K * (El) + (El) . (72) 

Unfortunately, the sudden-oscillation and sudden-decay approximations both become in- 
adequate in this regime, so that Eq. (72) will not in fact be correct. Departures in the 
sudden-oscillation can be calculated once the potential of the oscillating is known but no 
relevant calculation has been published so far. Even in the absence of such a calculation 
though, it follows just from continuity that that there will be a regime in which s c d m is within 
the present observational bound, while being big enough to be observable in the future. 



V. CONCLUSION 

The observed primordial curvature perturbation may well be accompanied by isocurva- 
ture perturbations in CDM, baryonic matter or neutrino. In the curvaton scenario, such 
a perturbation is inevitably generated by the curvaton perturbation, unless the curvaton 
decays before the relevant quantity (CDM, baryon number or lepton number) is created. 
This 'residual' isocurvature perturbation is fully correlated with the curvature perturbation. 

We presented a general expression, Eq. (46), which allows the residual isocurvature per- 
turbation to be evaluated once the mechanism for generating the relevant quantity is speci- 
fied. We then evaluated the residual CDM isocurvature perturbation for all the commonly 
considered candidates in the sudden decay approximation, and compared it with the ob- 
servational upper bound [31]. We first recalled previously-known results [17]. The CDM 
cannot be created right after inflation in the curvaton scenario, as in the case of CDM in 
the form of super- heavy weakly interacting massive particles (Wimpzillas), created from the 
vacuum energy [54,55]. But if the curvaton dominates the energy density before it decays, 
then the CDM can be created from that decay. 
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We then went on to consider new cases, using Eq. (47). After a fairly complicated 
calculation, we derive Eq. (63) for the case of CDM falling out of thermal equilibrium 
(weakly interacting massive particle or WIMP CDM). It is compatible with observation 
only if the CDM is created when the curvaton density has a rather specific value. Modifying 
the previous calculation slightly, we arrive at the case of axion CDM given by Eq. (69), 
which we find is too big in all cases. These results mean that axions, and probably also 
WIMPS, have to be created after curvaton decay. Such a requirement is not however a 
strong constraint on the curvaton scenario, since WIMPS and axions are both typically 
created rather late (at temperatures respectively of order 10 GeV and 1 GeV). 

Finally, we consider the case of CDM created from the oscillation of a scalar field with 
fixed mass. As the created particle would need to be massive and very weakly interacting to 
avoid subsequent thermalisation, we are again dealing with a Wimpzilla. We derive Eq. (70) 
for the residual isocurvature perturbation in this case. It again gives s c d m = _ 3 in the limit 
where the curvaton density when the oscillation starts is much smaller than its final value. 
In contrast with the previous cases though, it gives s c d m = in the opposite limit where 
curvaton domination is complete when the oscillation starts. There must therefore be an 
intermediate regime where s c d m in this case is of potentially observable magnitude. 

We can summarise the situation as follows. The residual CDM isocurvature perturbation 
is absent if the CDM is created after curvaton decay. Otherwise it is present, and it is 
typically either too big to be compatible with observation, or completely negligible. There 
are cases though where it can be of observable magnitude. In particular, this can occur if 
the CDM comes from the decay of the curvaton, or from the oscillation of a scalar field of 
fixed mass. In both of these cases, the residual isocurvature perturbation vanishes in the 
limit that the curvaton completely dominates the energy density at the time of the CDM 
creation, but isocurvature perturbations can exist at a detectable level if the domination is 
not complete. 

ACKNOWLEDGMENTS 

We are grateful to Y. C. Kim for a comment on our treatment of the WIMP 
case. This work was supported in part by PPARC grants PPA/G/S/1999/00138, 
PPA/G/S/2000/00115 and PGA/G/O/2000/00466. DW is supported by the Royal Society. 



18 



REFERENCES 



[1] E. W. Kolb and M. S. Turner, "The Early Universe," Redwood City, USA: Addison- 
Wesley (1990) 547 p. (Frontiers in physics, 69). 

[2] A. R. Liddle and D. H. Lyth, "Cosmological inflation and large-scale structure," Cam- 
bridge, UK: Univ. Pr. (2000) 400 p. 

[3] D. N. Spergel et al, arXiv:astro-ph/0302209. 

[4] H. V. Peiris et al, arXiv:astro-ph/0302225. 

[5] E. Komatsu et al, arXiv:astro-ph/0302223. 

[6] J. Maldacena, JHEP 0305, 013 (2003) [arXiv:astro-ph/0210603]; V. Acquaviva, N. Bar- 

tolo, S. Matarrese and A. Riotto, arXiv:astro-ph/0209156. 
[7] D. S. Salopek, Phys. Rev. D 45 (1992) 1139; N. Bartolo, S. Matarrese and A. Riotto, 
Phys. Rev. D 65, 103505 (2002) [arXiv:hep-ph/0112261]; F. Bernardeau and J. P. Uzan, 
Phys. Rev. D 67, 121301 (2003) [arXiv:astro-ph/0209330]. 
[8] D. H. Lyth and A. Riotto, Phys. Rept. 314, 1 (1999) [arXiv:hep-ph/9807278]. 
[9] D. H. Lyth and D. Wands, Phys. Lett. B 524, 5 (2002) [arXiv:hep-ph/01 10002]. 
[10] S. Mollerach, Phys. Rev. D 42, 313 (1990). 
[11] A. D. Linde and V. Mukhanov, Phys. Rev. D 56, 535 (1997). 
[12] G. Dvali, A. Gruzinov and M. Zaldarriaga, arXiv:astro-ph/0303591. 
[13] L. Kofman, arXiv:astro-ph/0303614. 

[14] K. Enqvist, A. Mazumdar and M. Postma, arXiv:astro-ph/0304187. 
[15] K. Dimopoulos and D. H. Lyth, arXiv:hep-ph/0209180. 

[16] D. Wands, N. Bartolo, S. Matarrese and A. Riotto, Phys. Rev. D 66, 043520 (2002) 

[arXiv:astro-ph/0205253]. 
[17] D. H. Lyth, C. Ungarelli and D. Wands, Phys. Rev. D 67, 023503 (2003) [arXiv:astro- 

ph/0208055]. 

[18] T. Moroi and T. Takahashi, Phys. Lett. B 522, 215 (2001). 
[19] N. Bartolo and A. R. Liddle, Phys. Rev. D 65, 121301 (2002). 

[20] T. Moroi and T. Takahashi, Phys. Rev. D 66, 063501 (2002) [arXiv:hep-ph/0206026]. 
[21] M. Fujii and T. Yanagida, Phys. Rev. D 66, 123515 (2002) [arXiv:hep-ph/0207339]. 
[22] A. Hebecker, J. March-Russell and T. Yanagida, Phys. Lett. B 552, 229 (2003) 

[arXiv:hep-ph/0208249]. 
[23] R. Hofmann, arXiv:hep-ph/0208267. 

[24] M. Bastero-Gil, V. Di Clemente and S. F. King, Phys. Rev. D 67, 103516 (2003) 

[arXiv:hep-ph/0211011]. 
[25] M. Bastero-Gil, V. Di Clemente and S. F. King, Phys. Rev. D 67, 083504 (2003) 

[arXiv:hep-ph/0211012]. 
[26] T. Moroi and H. Murayama, Phys. Lett. B 553, 126 (2003) [arXiv:hep-ph/0211019]. 
[27] K. Enqvist, S. Kasuya and A. Mazumdar, Phys. Rev. Lett. 90, 091302 (2003) [arXiv:hep- 

ph/0211147]. 

[28] K. A. Malik, D. Wands and C. Ungarelli, Phys. Rev. D 67, 063516 (2003) [arXiv:astro- 
ph/0211602]. 

[29] M. Postma, Phys. Rev. D 67, 063518 (2003) [arXiv:hep-ph/02 12005]. 

[30] B. Feng and M. Li, arXiv:hep-ph/0212213. 

[31] C. Gordon and A. Lewis, arXiv:astro-ph/0212248. 

[32] K. Dimopoulos, arXiv:astro-ph/0212264. 



19 



M. Giovannini, arXiv:hep-ph/0301264. 

A. R. Liddle and L. A. Urena-Lopez, arXiv:astro-ph/0302054. 
J. McDonald, arXiv:hep-ph/0302222. 

K. Dimopoulos, G. Lazarides, D. Lyth and R. Ruiz de Austri, JHEP 0305, 057 (2003) 
[arXiv:hep-ph/0303154]. 

K. Enqvist, A. Jokinen, S. Kasuya and A. Mazumdar, arXiv:hep-ph/0303165. 

K. Dimopoulos, D. H. Lyth, A. Notari and A. Riotto, arXiv:hep-ph/0304050. 

M. Endo, M. Kawasaki and T. Moroi, arXiv:hep-ph/0304126. 

M. Postma and A. Mazumdar, condensate," arXiv:hep-ph/0304246. 

S. Kasuya, M. Kawasaki and F. Takahashi, arXiv:hep-ph/0305134. 

D. H. Lyth & D. Wands, arXiv:astro-ph/0306498. 

K. Dimopoulos, G. Lazarides, D. Lyth and R. Ruiz de Austri, in preparation. 

D. Wands, K. A. Malik, D. H. Lyth and A. R. Liddle, Phys. Rev. D 62, 043527 (2000) 

[arXiv:astro-ph/0003278]. 

M. Bucher, K. Moodley and N. Turok, Phys. Rev. D 62, 083508 (2000) [arXiv:astro- 
ph/9904231]. 

R. Trotta, A. Riazuelo and R. Durrer, Phys. Rev. Lett. 87, 231301 (2001) [arXiv:astro- 

ph/0104017]; Phys. Rev. D 67, 063520 (2003) [arXiv:astro-ph/0211600]. 

L. Amendola, C. Gordon, D. Wands and M. Sasaki, Phys. Rev. Lett. 88, 211302 (2002) 

[arXiv:astro-ph/0107089]. 

J. Valiviita and V. Muhonen, arXiv:astro-ph/0304175. 

P. Crotty, J. Garcia-Bellido, J. Lesgourgues and A. Riazuelo, arXiv:astro-ph/0306286. 
D. Langlois, Phys. Rev. D 59, 123512 (1999). 

A. R. Liddle, D. H. Lyth, K. A. Malik and D. Wands, Phys. Rev. D 61, 103509 (2000) 
[arXiv:hep-ph/99 12473]. 

M. Sasaki and E. D. Stewart, Prog. Theor. Phys. 95, 71 (1996) [arXiv:astro-ph/9507001]. 
T. Damour and A. Vilenkin, Phys. Rev. D 53, 2981 (1996) [arXiv:hep-th/9503149]. 
D. H. Lyth, D. Roberts and M. Smith, Phys. Rev. D 57, 7120 (1998) [arXiv:hep- 
ph/9609441]. 

D. J. Chung, E. W. Kolb and A. Riotto, Phys. Rev. Lett. 81, 4048 (1998) [arXiv:hep- 

ph/9805473]; D. J. Chung, E. W. Kolb and A. Riotto, Phys. Rev. D 60, 063504 (1999) 

[arXiv:hep-ph/9809453]. 

D. H. Lyth, Phys. Rev. D 45 (1992) 3394. 

M. S. Turner, Phys. Rev. Lett. 59, 2489 (1987) [Erratum-ibid. 60, 1101 (1988)]. 
T. Moroi and L. Randall, Nucl. Phys. B 570, 455 (2000) [arXiv:hep-ph/9906527]. 



APPENDIX A: GAUGE-INVARIANT CALCULATION 

In this section we reproduce the calculation of the CDM perturbation given in section 
IV in terms of manifestly gauge invariant perturbations. 

To determine the CDM isocurvature perturbation we need to determine the gauge- 
invariant perturbation in the the CDM number density on spatially flat slices, Eq. (12), 
(or equivalently the curvature perturbation on uniform number-density slices) 
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~ _ 1 Sn cdm 

Ccdm = W- 

<J ^cdm 



(Al) 



where Sn cdm and ip are the gauge-dependent perturbed number density and curvature per- 
turbation respectively. 

In the sudden freeze- out approximation, we assume that the CDM remains in equilibrium 
abundance, with number density given by Eq. (54), until the local interaction rate per 
expansion time falls below the critical value 



r 

H 



K . 



(A2) 



Thereafter the CDM number is conserved and hence the curvature perturbation Ccdm in 
Eq. (Al) remains constant on large scales after freeze-out. The physical condition (A2) 
picks out a physical hypersurface £/ which marks the transition between equilibrium and 
freeze-out. Thus we can determine Ccdm in Eq. (Al) after freeze-out from the curvature 
perturbation, ipf, and perturbed number density, Snf, on this hypersurface. 

A uniform temperature slice at this time has no perturbation in the equilibrium 
number density, Sri? = 0, but may have a curvature perturbation ip? = — Cr- ( in the 
curvaton scenario we have Cr <C ( a before the curvaton decays.) We obtain ipf and Srif 
via a gauge-shift from the hypersurface £t to the hypersurface £/. This corresponds to a 
coordinate-shift 



We then have 



St 



fT 



8{T/H) T 



(T/H) - (HT/H 2 ) 



(A3) 



and hence 



ipf — ipT 4" HStfT , 
5rif = 5riT — nStfT , 



~ 1 (m 3\ 

Ccdm = Cr + - y— - -J H^dtfT , 



(A4) 
(A5) 



(A6) 



where we have assumed H~ X T /T — — 1. 

Assuming T/H oc T a /p 1 ^ 2 (where a > 2 for CDM to fall out of equilibrium with de- 
creasing temperature) we have from Eq. (A3) 



St 



4- a, 



fT 



\ Spr 

2{2-a)-n a )^ p* ' 



(A7) 



and using 



• Cr C — A n (Cr Co") j 



(A8) 



we then have 
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Finally we then obtain from Eq. (A6) and (A9) 




m 3\ Vt, 



(C - Cr) • 



(A10) 



T, 2/2(a-2) + ft, 



If we have adiabatic perturbations (( r = before freeze-out, we must have ( cdm = ( r 
from Eq. (A10), and hence s C( jm — after freeze-out. 

In the curvaton scenario when the CDM decouples before the curvaton decays, we have 
Cj ^> Cr and the primordial curvature perturbation is given by ( = r( a . The residual CDM 
isocurvature perturbation is then given by Eq. (63) of section III. 
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